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Abstract 

The dynamics of the eigenvalues (semimartingales) of a Levy process X with values in Hermitian 
matrices is described in terms of ltd stochastic differential equations with jumps. This generalizes the 
well known Dyson-Brownian motion. The simultaneity of the jumps of the eigenvalues of X is also studied. 
If X has a jump at time t two different situations are considered, depending on the commutativity of 
X{t) and X{t—). In the commutative case all the eigenvalues jump at time t only when the jump of X 
is of full rank. In the noncommutative case, X jumps at time t if and only if all the eigenvalues jump at 
that time when the jump of X is of rank one. 

Key words/ Dyson-Brownian motion, infinitely divisible random matrix, Bercovici-Pata bijection, 
matrix semimartingale, simultaneous jumps, non-colliding process, rank one perturbation, stochastic 
differential equation with jumps. 
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1 Introduction 


This paper is on a topic that meets a number of subjects where Ole Barndorff-Nielsen has contributed. 
One is Ole’s interests in the study of infinitely divisible random matrix models and their associated matrix¬ 
valued processes, like subordinators, Levy and Ornstein-Uhlenbeck processes Another is the 

connection between classical and free infinite divisibility, where Ole’s papers with Thorbjprnsen m-m 
drove the study of the so-called Upsilon transformations of classical infinitely divisible distributions m, i, 

El- 

Matrix Levy processes with jumps of rank one arise as a covariation process of two C'’*—valued Levy 
processes |23] and they provide natural examples of matrix subordinators and matrix Levy processes of 
bounded variation, such as those considered in [3] and [S]. 

In the context of free probability, Bercovici and Pata [14] introduced a bijection from the set of classical 
infinitely divisible distributions on R to the set of free infinitely divisible distributions on R. This bijection 
was described in terms of ensembles of Hermitian random matrices : d > 1} by Benaych-Georges [T3] 
and Cabanal-Duvillard For fixed d > 1, the associated d x d matrix distribution of is invariant 
under unitary conjugations and it is (matrix) infinitely divisible. The associated Hermitian Levy process 
{Xd{t) : t>0} has been considered in [23] and it has the property that, in the non pure Gaussian case, the 
jumps AXd{t) = Xd[t) — Xd(t—) are d x d matrices of rank one, see also [52], [55] . 

The simplest example of this connection is the well-known theorem of Wigner that gives the semi¬ 
circle distribution (free infinitely divisible) as the asymptotic spectral distribution of {Xd : d > 1} in the 
case of the Gaussian Unitary Ensemble (GUE), see [T]. In this case, for each fixed d > 1, the associated 
Levy process of Xd is the d x d Hermitian Brownian motion {B{t) : t > 0} = {{Bjk{t)) : t > 0} where 
(Hjj(t))^^j^, (ReHjfc(t))j<fc, (Im Bjk{t))j<^k is a set of d^ independent one-dimensional Brownian motions 
with parameter |(1 -I- d^fe). 

Still in the GUE case, let A = {(Ai(<), X 2 {t),Xd{t)) : t > 0} be the d-dimensional stochastic process of 
eigenvalues of H. In a pioneering work, Dyson |24] showed that if the eigenvalues start at different positions 
(Ai(0) > A2(0) > ... > Aci(O) a.s.), then they never meet at any time (Ai(t) > X 2 {t) > ... > Xd{t) V t > 0 
a.s.) and furthermore they form a diffusion process (and a semimartingale) satisfying the ltd Stochastic 
Differential Equation (SDE) 

dA,(f) = dWi{t) + Y^ A-(f) -V(f) ’ ^^0,l<z<d, (1) 

where Wi, ...,Wd are independent one-dimensional standard Brownian motions (see [T] Sec. 4.3.1]). The 
stochastic process {X{t) : t > 0} is called the Dyson non-colliding Brownian motion corresponding to the 
GUE, or, in short, Dyson-Brownian motion. The study of the eigenvalue process of other continuous matrix¬ 
valued processes has received considerable attention, see m -[19]. [25]-[28]. The eigenvalues of all these 
processes are strongly dependent and do not collide at any time almost surely. 

The aim of this paper is to understand the behavior on time of the eigenvalue process A = {(Ai(t), A 2 (f),... 
,Xd{t)) : t > 0} oi a. d X d Hermitian Levy process X = {A(t) : f > 0}. The first goal is to give conditions 
for the simultaneity of the jumps of the eigenvalues. In particular, when the jump of X is of rank one and 
X{t) and X(t—) do not commute, we show that a single eigenvalue jumps at time t if and only if all the 
eigenvalues jump at that time. This fact which is due to the strong dependence between the eigenvalue 
processes contrasts with the case of the no simultaneity of the jumps at any time of d independent real Levy 
processes; in Example [6] we give a Hermitian Levy process with such eigenvalue process. We observe that in 
general, the eigenvalue process is not a Levy process (as in the Dyson-Brownian case) but a semimartingale 
given in terms of the Hermitian Levy process X. 

The second goal is to describe the dynamics of the eigenvalue process, analogously to O- This leads to 
first considering the appropriate Ito’s formula for matrix-valued semimartingales with jumps, more general 
than those of bounded variation considered by Barndorff-Nielsen and Stelzer [8] . 

Our main results and the structure of the paper are as follows. In Section 2 we gather several results on 
Hermitian Levy processes, including versions of Ito’s formula which are useful for studying the corresponding 
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eigenvalue processes. We also point out conditions under which the spectrum of X will be simple for every 
t > 0 and conditions for the differentiability of the eigenvalues. In Section 3 we consider the simultaneous 
jumps of the eigenvalue (semimartingale) process of a d x d Hermitian Levy process X. When AX{t) = 
X{t) — X{t—) is not zero and X{t) and X{t—) commute, we show that the eigenvalues jump simultaneously 
if and only if the matrix jump AX{t) is of rank d. On the other hand, if AX{t) 0 is of rank one and X{t) 
and X{t—) do not commute, then the process X jumps at time t > 0 if and only if all the components of A 
jump at that time. This leads to analyzing the relation between the eigenvalues of X{t) and X{t—) when the 
jump AX{t) is a rank one matrix. A key step in this direction is the recent paper |31] on the spectrum of a 
rank one perturbation of an unstructured matrix. Extensions to Hermitian semimartingales are also pointed 
out. Finally, Section 4 considers the dynamics of the eigenvalues of a Hermitian Levy process, extension to 
the Dyson-Brownian motion, showing that a repulsion force appears in the bounded variation part of the 
eigenvalues only when there is a Gaussian component of X. One of the main problems to overcome is the 
fact that even when X{t) is Hermitian and has a simple spectrum in an open subset almost surely, X{t—) 
may not. 

We shall use either X{t) or Xt to denote a stochastic process in a convenient way when the dimension 
or the entries of a matrix have also to be specified. 


2 Preliminaries on Hermitian Levy processes 


In this section we gather several results on matrix-valued semimartingales and general Hermitian Levy 
processes. 

Let Md = Md (C) denote the linear space of d x d matrices with complex entries with scalar product 
{A,B) = tT{B*A) and the Frobenius norm ||A|| = [tr(A*A)]^^^ where tr denotes the (non-normalized) 
trace. The set of Hermitian random matrices in Md is denoted by Hd, H)) = ]HId\{0} and H;) is the set of 
rank one matrices in Hd. 

We will use the identification of a d x d Hermitian matrix X = (xij + iyij)i<ij<d with an element in , 
that is. 


X — {^ij “f iyij )l<z,j<d t-)- X dzd)2<i<j<d^ 


( 2 ) 


Thus, the set Hd is identified with and one has that for X,Y € Hd, tr (XY) = Y^X is real. 

An Md-valued process X = {{xij){t) : t > 0} is a matrix semimartingale if Xij{t) is a complex semi¬ 
martingale for each i,j = l,...,d. Let X = {{xij){t) : t > 0} G Hd and Y = {{yij){t) : t > 0} G Hd be 
semimartingales. Similar to the case of matrices with real entries in [8], the matrix covariation of X and Y 
is defined as the Hd-valued process [A, T] := {[A, T] (t) : t > 0} with entries 


d 

W = [x^k,yk]] {t), 




(3) 


where [xik^ykj] (t) is the covariation of the C-valued semimartingales {xik{t) : t > 0} and {ykj{^) ■ t ^ 0}, 
see pp 83]. One has a decomposition into a continuous part and a pure jump part by 

[A, A] it) = [A^ y^j it) + AA« Ay;, (4) 

s<t 


where [A; y% (t) := 


x^k^yko 


(t). We recall that for any semimartingale a;, the process is the a.s. 


unique continuous local martingale m such that [x — m] is purely discontinuous. 

Matrix-valued Levy processes are examples of matrix-valued semimartingales. In particular, so are the 
Hermitian Levy processes. 
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2.1 The Levy—Khintchine representation and Levy—Ito decomposition 

A d X d matrix Hermitian Levy process X = {X{t) : t > 0} is characterized by the Levy-Khintchine repre¬ 
sentation of its Fourier transform = exp(t'0(0)) with Laplace exponent 

^(0) = itr(0vl/ ) - itr (0.40) + f ( - 1 - ) ^^(dO, O e 

^ -'Hrf y i-f iicii / 

where A : ^ is a positive symmetric linear operator (i.e., tr ($.4$) > 0 for $ e and tr (02.40i) = 

tr (01.402) for 01,02 G Mji^), is a measure on Hd (the Levy measure) satisfying i^({0}) = 0 and /jj_j(l A 

||a;||^)^(da:) < oo, and G Hd. The triplet uniquely determines the distribution of X. 

Interesting matrix Gaussian distributions for some examples of A are the following (here '1 = 0 and 
iy = 0): 

(i) A = aAd, > 0, corresponds to the GUE case of parameter cr^; 

(ii) More generally, .40 = Ei 0 E 2 , for Ei, E 2 nonnegative definite matrices in H^, corresponds to the 
matrix Gaussian distribution with Kronecker covariance Ei 0 E 2 ; 

(iii) .40 = tr(0)cr^ld, cr^ > 0, is the covariance operator of the Gaussian random matrix gld where g is 
a one-dimensional random variable with zero mean and variance cr^. 

Any Hd-valued Levy process X = {A(<) : t > 0} with triplet (A, :^,'!) is a semimartingale with the 
Levy-Ito decomposition 

X{t) = t'i + + [ f yJx{ds,dy)+ [ [ 2 / 4 x(ds, dy), t > 0 , (5) 

Jl0,t]J\\y\\<l dlO,t]J\\v\\>l 

where : t > 0} is a Hd-valued Brownian motion with covariance A, i.e., it is a Levy process with 

continuous sample paths (a.s.) and each Bj\^{t) is centered Gaussian with 

E {tr(0ii?^(s))tr (02l3.4(t)) } = min(s, t)tr (0iA02) for each 0i, 02 G Hj;; 

Jx(', •) is a Poisson random measure of jumps on [0, 00 ) x Hj) with intensity measure Leb ( 8 > v, independent 
of {Bj^{t) : t > 0 } and the compensated measure is given by 

Jx{dt,dy) = Jx{dt,dy) -dtv{dy). 

For a systematic study of Levy processes with values in Banach spaces, see [2]. 

2.2 Hermitian Levy processes with simple spectrum 

Hermitian Levy processes with simple spectrum (all eigenvalues distinct) are of special importance in this 
paper. For each t > 0, if X{t) has an absolutely continuous distribution then X{t) has a simple spectrum 
almost surely, since the set of Hermitian matrices whose spectrum is not simple has zero Lebesgue measure 
[54] . see also [Tj. The absolute continuity gives also that X{t) is a nonsingular random matrix since for a 
Hermitian matrix A, det(A) is a real polynomial in variables and {A G Hd : det(A) = 0} has zero Lebesgue 
measure in . 

If there is a Gaussian component, i.e., A 7 ^ 0 , then X{t) has an absolutely continuous distribution for each 
< > 0. Next we point out a sufficient condition on the Levy measure of a Hermitian Levy process without 
Gaussian component to have an absolutely continuous distribution and hence a simple spectrum for each 
fixed time. 

Let Sd be the unit sphere of H^. A polar decomposition of the Levy measure 1 / is a family of pairs (tt, p^) 
where tt, the spherical component, is a finite measure on Sd with 7 r(S£i) > 0 and the radial component, 
is a measure on (0, 00 ) for each ^ G Sd with P{((0, 00 )) > 0 such that 

,^{E) = [ nidO [ lEiuOpd^u), E G B(M°d). 

JSd (0,oo) 
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We say that v satisfies the condition D if, for each ^ G Sd, there is a nonnegative function on (0, oo) 
such that (B) = Jg (u) du for any B G B(M\ {0}) and the following divergence condition holds 


/ g^ {u) du = oo n-a.e. 

J (0,00) 

As an example, recall that a random matrix M is called self-decomposable (and hence infinitely divisible) 
if for any b G (0,1) there exists Mf, independent of M such that M and bM -|- Mb have the same matrix 
distribution. An infinitely divisible random matrix M in with triplet is self-decomposable if 

and only if 

iy{E) = [ 7T{dO f lEio'^idu), EgB{E°), 

JSd d ( 0 ,oo) ^ 

with a finite measure tt on and a nonnegative decreasing function on (0,oo) for each ^ G Sd- In this 
case u satisfies the condition D for g^{u) = k^{u)/u where we can choose the measure tt to be zero on the 
set e Sd : fc^(O-l-) = 0}. 

The following proposition for Hermitian Levy processes follows from the corresponding real vector case, 
Theorems 27.10 and 27.13 in [51], by using the identification (Hj) of Hd with K'^^. We recall that a Levy 
process is nondegenerate if its distribution for any fixed time is nondegenerate. 

Proposition 1 (a) If {X{t) : t > 0} is a nondegenerate Levy process in Hd without Gaussian component 
satisfying condition D, then Xft) has absolutely continuous distribution for each t > 0. 

(b) If {X{t) : t > 0} is a nondegenerate self-decomposable Levy process in Hd without Gaussian compo¬ 
nent, then X{t) has absolutely continuous distribution for each t > 0. 


2.3 Smoothness of the spectrum of Hermitian matrices 

Next we consider several facts about eigenvalues as functions of the entries of a Hermitian matrix, similar 
to |28j for symmetric matrices, that will be useful in the sequel. 

Following [I], let be the set of unitary matrices U = (uij ), with uu > 0 for all i, utj ^ 0 for i, j, and 
all minors of Lf having non-zero determinants. Let us denote by H^^ the set of matrices X in Hd such that 


X = UDU*, 


( 6 ) 


where D is a diagonal matrix with entries \i = Da such that Ai > A 2 > • • • > Ad, 1/ G An element of 

H^*^ is called a very good matrix. The set of very good matrices H^^ can be identified with an open subset 
of . It is known that the complement of H^*^ has zero Lebesgue measure. If T : —)• is the 

mapping defined by 


rr/TT\ — I ^12 Uld U23 U2d M(d-l)d 

1 \U ) - I 1 ’ • ' 1 1 5 • • • ? 1 • ’ • 1 

uii un U22 U22 U(d-i)(d-i) 


then Lemma 2.5.6 in [I] shows that T is bijective, smooth, and the complement of is a closed subset of 
zero Lebesgue measure. 

Let Sd be the open set 


Sd = {(Ai, A2,..., Ad) G : Ai > A2 > • • • > Ad} 


(7) 


and let D\ be the diagonal matrix Du = Ai for any A = (Ai, A 2 ,..., Ad) G Sd- Then the map T : Sd x 
T (U^'^) Elddefined by T (A, z) = T~^{z)*DxT~^{z) is a smooth bijection and the inverse mapping of 
T satisfies _ _ 

T-i(X)= (a(X),C/(X)), 
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where U{X) = T{U) for U € satisfying ([5]). As a consequence of these results, there exist an open 
subset G C such that \G has zero Lebesgue measure and a smooth function 

$ : ( 8 ) 

such that A (A) = (A) for A e G. 


2.4 The Ito Formula for Hermitian Levy processes 


We shall first consider a convenient notation for the Frechet derivative of a function / from Hd to R 
(see m X.4]), using the identification ([2]) of an element A = {xij + iyij)i<i,j<d in Eld as the element 

^xii, ...,XddT in R . 

The derivative of / at the point u is denoted by Df{u); it can be thought as a linear operator from R'^ 
to R. Its representation as a Hermitian matrix is 


Df 


dxi 


^y^J 


l<i<j<d 


(9) 


Let L ^R'^^jRj denote the space of linear operators from R'^^ to R. If the map Df : R'^^ —>• L ^R'^^jRj 

defined by u —>■ Df(u) is differentiable at a point u, we say that / is twice differentiable at u. We denote 
the derivative of the map Df at the point u by D‘^f{u). Recall that D^f(u) is a bilinear operator from 
Rd X R'^ to R which is symmetric, that is, D'^f{u){x,y) = D^f{u){y,x). We use the following notation for 
the first and second partial derivative operators, where Zij = Xij or Zij = ytj, with ya = 0, i = 1, ...,d 


T. 





l<i<d 


dXii 

l<i<j<d 


+ 


dyi^ 


( 10 ) 


E 


dzij Zkh 


= E 

l<i<j<d 




dXiiXjj 


+ E E 

l<k<.h<d 


+ E E ( 

l<k<.hKd 


d’^ 

dxijXkh 


dyijykh' 


dxijykh 


) 


( 11 ) 


For a linear operator A : Hd —^ H^, we denote its dP x df matrix representation by Al = {Aij^kh)- 

Following Barndorff-Nielsen and Stelzer [8], we say that for each open subset G C the process 
A = {At : t > 0} is locally bounded within G if there exist a sequence of stopping times (T„)^>^ increasing 
to infinity almost surely and a sequence of compact convex subsets Dn C G with C -Dn+i such that 
At G Dn for all 0 < t < T„. Observe that if the process A is locally bounded within G, it cannot get 
arbitrarily close to the boundary of G in finite time, hence At and At_ are in Dn at all times t > 0. 

The following result can be proved similarly to Proposition 3.10.10 in Bichteler [16], who considered the 
multivariate case. The bounded variation multivariate case can be seen in [^. 


Theorem 2 (The Ito formula for matrix semimartingales) Let (At : t > 0} 6e a matrix semimartin¬ 
gale, G C Eld an open subset, and / : G ^ R continuously differentiable. If {At : t > 0} is locally bounded 
within G, then At_ belongs to G for all t > 0, the process f (At) is a semimartingale, and 

f (At) = / (Ao) + tr f /* Df (A«_) dA,) + 1 ^ i^s) d[zt^,zl^]s 

\Jo+ J ^ Jo+ 


+ f [f (A,) - / (A,_) - tr [Df (A,_) y)] Jx {ds, dy ). 

d(0,t]xH« 
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The following two propositions for Hermitian Levy processes are special cases of the above result. The 
first one is the analog of the multivariate case in Prop. 8.15]. 

Proposition 3 (The Ito formula for matrix Levy processes) Let {Xt '■ t > 0} be a Hermitian Levy 
process with triplet G C Hj; an open subset, and / : G —>■ R twice continuously differentiable. 

If {Xt : < > 0} is locally bounded within G, then Xt- belongs to G for all t > 0, the process / {Xt) is a 
semimartingale, and 

f (Xt) = f (Xo) + tr ( r Df (X,_) dXs) f (X,) ds (12) 

^ {i,j){k,h) OZzjZkh 


+ f [f (Xs- +y)- f (Xs-) - tr {Df {Xs-) y)] Jx {ds, dy). 

J(0,i]xHO 

The next martingale-drift decomposition of functions of a Hermitian Levy proces X can be obtained by 
replacing the Levy-Ito decomposition ([5]) for X into the stochastic integral part of (fT^ . 


Proposition 4 Let {Xt : t > 0} be a Hermitian Levy process with triplet {A, G C an open subset, 

and / : G —^ R twice continuously differentiable. If {Xt : t > 0} is locally bounded within G, then Xt- 
belongs to G for all t > 0, and the process f {Xt) can be decomposed as f {Xt) = Mt + Vt where Mt is a 
martingale and Vt is a bounded variation process such that 

Mt = f (Xo) + tr f / Df {Xs-) BA{ds)) + f [f (X,_ + y) - f (X,_)] Jx {ds,dy) 

\Jo / d(o,t]xH“ 

and ^ ^ 

Xij^khjr^ - f{Xs)ds + tr([ 4'H/(Xs)ds 

2 -^0 

+ [ [f {Xs- +y)- f {Xs-) - tr {Df {Xs-)y) l{||y||<i}] dsv{dy). 

d(0,t]xH“ 


Remark 5 If a Hermitian Levy process X is invariant under unitary conjugation with triplet (^cr^lj^ 2 , v, fAd), 
where v has the polar decomposition (tt, p), and tt is the uniform measure on Sd, then the ltd formula for X 
in Proposition becomes in 


Mt = f (Xo) + tr ^ Df (X«_) Bj^{ds)^ 

[/ (Xs_ +ru) - f (Xs_)] {Jx {ds, du, dr) - ds'K{du)p{dr)) 


' (0,t] xSd X (0,oo) 


and 


Vt = 


E 

{i,j){k,h) 


52 


dzij Zkh 


f {Xs) ds + V'tr Df {Xs) ds 


+ [ [f {Xs-+ru) - f {Xs-) - rti {Df {Xs-)u)l[,.<iy] dsTT{du)p{dr). 

J (0,i] xSd X (0,oo) 
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3 Simultaneous jumps of a Hermitian Levy process and its eigen¬ 
values 

Let : t > 0} he a. d X d Hermitian Levy process and let A(t) = (Ai(t),Ad(t)) be the vector of 

eigenvalues of X{t) where Ai(t) > A 2 (t) > ••■ > Ad(t) for each t > 0. If A jumps at time t, it is clear 
that AX{t) 7 ^ 0, since the d-tuple of ordered eigenvalues is a continuous function on the space of Hermitian 
matrices. This follows from the Hoffman-Wielandt inequality (see [TJ Rem. 2.1.20]) where for each t >0 

||AA(t)||2<||AX(t)||, (13) 

where || II 2 denotes the Euclidean norm in 

A natural question is whether the jump of the process X at time t implies the jump of the eigenvalues 
at that time. Our goal in this section is to provide conditions for the simultaneous jumps of the eigenvalue 
process A whenever X jumps. 

We analyze the non zero jumps of X in two situations, when X{t) and X(t—) commute and when they 
do not commute. In the commutative case we find that exactly k eigenvalues jump when X has a jump of 
rank k at time t, and that all the eigenvalues would jump simultaneously if and only if the jump of X at 
time t is of rank d. In the noncommutative case, we prove that all the eigenvalues jump simultaneously if 
the jump of X at time t is of rank one. 

We present below natural examples of both, commutative and noncommutative situations. 

Example 6 Let {(Ai(t),..., Ad(t)) : t > 0} be a d-dimensional Levy process where {Ai(t) : t > 0}, i = I, ...,d, 
are independent one-dimensional Levy processes. For an arbitrary non random unitary d x d matrix U, 
define the d x d Hermitian Levy process 

X{t) = Lfdiag{Xi{t ),..., Xd(t))U*, t >0. 

Since independent Levy processes cannot jump at the same time almost surely, X is an Hermitian Levy 
process whose jumps are of rank one almost surely and we observe that in that case X(t) and Xft—) are 
simultaneously diagonalizable. 

Example 7 Let us consider a compound Poisson d x d matrix process 

N{t) 

X(t) = UjU*, t>0, 

1=1 

where Uj,j > 1 are i.i.d. random vectors uniformly distributed on the unit sphere of C‘^ and N{t),t > 0, is 
a Poisson process independent of Uj,j > I. Since AXft) = niv(i)it^(t)l{AX(t)=i} the process X has jumps 
of rank one almost surely and Xlt) = Xft—) + 'UAr(t)U^(t) whenever there is a jump at time t. Hence we 
observe that Xlt) and X{t—) do not commute since and Xlf—) = uiu\ 

not commute. 

Example 8 Let us consider {X{t) : t > 0} a Hermitian Levy process in H])' the set of nonnegative definite 
d X d Hermitian matrices. Assume that X has jumps of rank one. Then by Theorem 4.1 in X is the 
quadratic variation of a -valued Levy process {Y{t) : t > 0}, that is 

X{t) = [V] (t), t > 0. 

For each t > 0, X{t) = ^ AW (AW)* and X{t—) = lim ^ AW (AW)* do not commute in general. In 

0<u<t ^F'*0<r<s 

particular, X{t) and X{t—) commute if for each u,r G (0, t] the term AW (AW)* AW (AW)* is Hermitian. 


Example 9 Let {Z(t) : t > 0} be a Hermitian Levy process in Htj of bounded variation with jumps of rank 
one. By Theorem 4.2 in \2S]I Z is the difference of the quadratic variations of two C‘^-valued Levy processes 
{X{t) :t>0} and {Y{t) :t>0}, 


Z{t) = [X] {t) - [F] (0, t > 0, 

where [X] {t),t > 0 and [Y] {t),t > 0 are independent Levy processes in It follows from Examplel^ that 
Z{t) = [X] (f) — [F] (i) and Z(t—) = [X] (<—) — [F] (t—) do not commute in general. In particular, Z{t) and 
Z{t—) commute if for each u,r € (0,t] the terms (AX„)* A (X^)* , AF^ (AF„)* A Yr (AF^)* and 
AXu (AXu)* A Yr (AFj.)* are Hermitian. 

3.1 Commutative case 

Assume that X(t) and X(t—) commute. Then there exists a random unitary matrix U = Ut such that 
X{t) = Udiag{Xilt),..., Xd{t))Lf* and X{t—) = Udiag{Xi{t—),..., Xd{t—))U*. Observe that the eigenvalues 
of AX{t) = X{t) — X{t—) are the difference between the eigenvalues of X{t) and the eigenvalues of X{t—). 

If AX{t) is of rank one, there is exactly one eigenvalue Xi such that AX{f) = Udiag{0,..., AXi{t),..., 0)1/* 
and all the other eigenvalues do not jump at that time t. Similarly if AX(t) is of rank k where 1 < fc < d 
there are exactly k eigenvalues that jump at time t and all other d — k eigenvalues remain without jumping 
at t. Observe that the only case when all eigenvalues jump simultaneously at t is when the jump AX{t) ^ 0 
is of rank d. 

3.2 Noncommutative case 

When X{t) and X{t—) do not commute the situation is completely different from the commutative case, 
since the eigenvalues of X{t) — X{t—) are not any more the difference between the eigenvalues of X{t) and 
the eigenvalues of X{t—). In this case we prove that if the jump AX{t) is of rank one, all eigenvalues jump 
simultaneously at this time t, as we will see in Theorem 1121 

We first consider some facts about the eigenvalues of rank one matrix perturbations. Let A he a d x d 
complex matrix. A rank one perturbation of A is a d x d complex matrix A + ruv^ where u,v G C‘^ and 
r G C. If S' is any set of polynomials in C[a;i, ...,Xd\, we define V{S) = {u e C''* : f{u) = 0 for all / G S}. 
If S is a finite set of polynomials /i,...,/„ we write F(/i,...,/„) for V{S). Recall that a subset W of 
is an algebraic set if IF = F(/i, ...,/„) for some /i, ...,/„. A non empty subset Q C is generic if the 
complement C‘^\Q is contained in an algebraic set IF which is not C''*. In such a case, C‘^\Q is nowhere 
dense and of 2d dimensional Lebesgue measure zero. 

When a Hermitian Levy process X has simple spectrum for each t > 0, then Xt and Xt- have no common 
eigenvalues when the jump at time t is of rank one. 

Proposition 10 Let {Xt : t > 0} be an absolutely continuous d x d Hermitian Levy process. Let X{t) = 
(Ai(t),..., Ad(t)) be the vector of eigenvalues of Xt for each t >0. If X has a rank one jump at time t > 0, 
then 

{Ai(W),..., Ad(Xt)} n {Ai(W_),..., Xd{Xt-)} = 0 a.s. 

Proof. Observe that Q = {z G : tr (zz*) ^ O} is a generic set since [z G : tr (zz*) = 0} = {0} is an 
algebraic set which is equal to F(/i,..., f 2 d) where fi{z) = zi for any z = (zi,..., Z 2 d)^ G and i = 1,..., 2d. 

Since X has a jump of rank one at time t, AXt = ruu* with r = rt G R\{0} and u = Ut G C''*\{0}. That 
is, we have the rank one perturbation Xt- = Xt + ruu*. Since Xt is absolutely continuous and Hermitian, 
the eigenvalues Ai(Xt), • ■ • , XdiXt) do not collide almost surely, hence the degree of the minimal polynomial 
of Xt is d. By Proposition 2.2 (hi) of [3T] applied to z = {u,u) G Q, there are exactly d eigenvalues of Xt- 
that are not eigenvalues of Xt. ■ 
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Remark 11 Proposition 2.2 in m also gives insight into the behavior of the non-Hermitian rank one 
jumps of a general complex matrix Levy process for which Xt no longer has a simple spectrum. In this case, 
AXt = ruv^ with r G C\{0} and u,v G the number mt < d of distinct eigenvalues of Xt might be 

random and there are exactly mt eigenvalues of Xt- that are not eigenvalues of Xt- 

We now present the main result of this paper, giving conditions for the simultaneous jumps of a Hermitian 
Levy process and its eigenvalue process. Moreover, if one eigenvalue component of the eigenvalue process 
jumps, then all the other eigenvalue components jump. 

Theorem 12 Let {X{t) : t > 0} be an absolutely continuous d x d Hermitian Levy process and let X{t) = 
(Ai(t),..., Xd{t)) be the vector of eigenvalues of X(t) where Xi(t) > X 2 {t) > • • • > Xd{t) for each t>0. 

If X has a jump of rank one at t > 0 and X{t) and X(t—) do not commute then AX(t) 7 ^ 0. Moreover 
AAi {t) 0 for some i if and only if AXj {f) 7 ^ 0 for all j 7 ^ i. 

Proof. If the Levy process X jumps at time t, then, by Proposition [TUI Xi{Xt) — Xj{Xt-) 7 ^ 0 for all 
i,j = 1 , 2 ,..., d and therefore AX{t) 7 ^ 0. 

Moreover, if A jumps at time t, it is clear from (fTUll that AX{f) 7 ^ 0, and then by the first part of the 
proof all the eigenvalues jump, that is AXj (t) 7 ^ 0 for all j ^ i. ■ 

Remark 13 a) We finally observe that Proposition 1101 and Theorem \1S\ are true also for a Hermitian 
semimartingale X with simple spectrum for each t>0. This is guaranteed whenever X{t) has an absolutely 
continuous distribution for each t > 0. 

b) However, we have written the results of this section in the framework of Hermitian Levy processes 
since the subject is one of Ole’s interests. 


4 The dynamics of the eigenvalues of Hermitian Levy processes 

4.1 Eigenvalue derivatives of Hermitian Levy processes 

Consider a Hermitian matrix X with decomposition X = UDU* given by ([U]) where U = (uij) is a unitary 
matrix, with ua > 0 for all i, Uij 7 ^ 0 for i,j, and all minors of U having non-zero determinants and D is a 
diagonal matrix with entries Xi = Du such that Ai > A 2 > ■ • ■ > Ad- 

Now if we assume that X = {xij + iyij)i<ij<d is a smooth function of a real parameter 9, the Hadamard 
variation formulae hold: for each m = 1 , 2 ,..., d, 


deXm = {U*deXUX 


doXra = {U*d^,XU)^^ + 2Y, 

j^m 


\{U*dgXU)^ 


Am X, 


In particular, for 9 = Xkh and ykh with k < h, we observe that 

dXr, 


dx 


kh 


— 2 R,6(zi/cm‘^hm ) “1“ 2|Lifc77j,| l{fc—/i}, 


dX,n 

dykh 


2 Im(zifcm‘a/im) ; 


Am 


= 2E 

j^m 


\ukmUhj UfimUkjl 

Xjn Xj 


■1 


{k^h} 


j^m 


\ukmy^kj I 


^k—h}, 
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9yL 


\'^k7n‘^hj Uhm‘^kj\ 

= ^ > -^^- k{k^h}. 


j^m 


Xm. Xq 


We point out that the spectral dynamics of Hermitian matrices is studied in [33], giving a Hadamard 
variation formula of any order for each eigenvalue in terms of a gradient. It is noted there that an eigenvalue 
depends smoothly on X whenever that eigenvalue is simple (even if other eigenvalues are not). 

Using the former derivatives of the eigenvalues and the smooth function given in (|8|) we have the following 
result for Hermitian Levy processes with simple spectrum for each time t > 0. 


Lemma 14 Let {Xt '■ t > 0} be a d x d Hermitian Levy process with absolutely continuous distribution for 
eacht > 0. Let X{t) = {Xi{t),Xd{t)) with Xi{t) > ••• > Xd{t) the eigenvalues of Xt for each t > 0. For each 
TO = 1, 2,..., d, there exists R which is C°° in an open subset G C , with Leb \G^ = 0, 

such that for each t > 0, Xm{t) = (Xt) and for k < h 


d<^>n 

dx 


kh 




(14) 


dykh 


— ^2 linl^UkmUhm')^{k^h} : 


and 


d<H 


= 2 


E \ukmUhj + UhmUkjl ^ n \ukmUkj\ ^ 

-- ifMM + 22^ X _w 1 


j^m 


dy. 


Xm X 

= 2E 

j^m 


I'^km’^hj '^hm'^kj 
Xm. Xq 


j^m 

|2 






(15) 

(16) 
(17) 


where for simplicity, the dependence on t is omitted in Xm = Xm{t) and Uki = Uki{t), which are random. 


4.2 The SDE for the process of eigenvalues 

Using Lemma[T4]and PropositionjS] we can describe the dynamics of the process of eigenvalues of a Hermitian 
Levy process. 


Theorem 15 Let {Xt : t > 0} be a Hermitian Levy process with triplet (.4, v, $) and locally bounded within 
the open set G in identified with the open set G in of Lemma \14\ Assume that Xt has an absolutely 
continuous distribution for each t > 0. Let X{t) = (Ai(t),..., Xd{t)) be the eigenvalues of Xt for each t > 0. 
Then, for each m = 1,... ,d the process Xm {Xt) is a semimartingale and 


Am (4t) — Am (4o) + tr r / ZlAm (4s_) + - f Aij kh-^ - Xm{Xs)ds 

\Jo J 2^0 dzqjZkh 


(18) 


+ 


'(0.t]xH° 


(i,j)(k,h) 

[Am {Xs- +y)- Xm {Xs-) - tr (DXm (X^.) ?/)] Jx {ds, dy) 


where a) 


DXm (Xs-) — 

9 $m(s-) 


9 $m(s-) , 9 $m(s-) 


dxkh 


+ 1 - 


dykh 


-1 


{k^h} 


l<k<h<d 


dxkh 


— 2Re(Ufem(s )uhm(s )) ^-{k^h} F 2 \ukm{s )| l{fc=?i}, 


(19) 
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h) 


d^m{s-) 

dykh 


2 Im {Uk m {s-)uh 771 (s-)) 


E 


32 ^ 

-^ij,kh 'o ] (-Eg) — 2 ^ ^ 

CfZijZkh 


E 


(^) (^) 


+ 4 L ^ 

l<fc</i<(i 


kh,kh ^ ^ 

j^m 


\Ukmis)uhjis)\ + |Uftm( 5 )Mfci(s)r 
(^) 


( 20 ) 


( 21 ) 


Proof. By Lemma [Til for each m = 1, 2,..., d there exists ffi. which is C°° in an open subset G 

of such that Am(i) = (^t) for each t > 0 , and the satisfy (fTH) - (fT71) . 

Now, taking / = in Proposition [31 where the open set G is identified with the open subset G, we 
have that Xm (Xt) is a semimartingale where formula (113]) becomes (fT3l) . The derivative matrix DXm (3fs-) 

in (dH), see ([9|), is given by DXm (Xg-) = ( ’ where the expressions of 

™ ^'^‘4 (ITSi) turn out to be in (IT^ and (l 20 l) . 

The second derivative term in (dH]) is given by J^ijM dz^z^^ Xm (Xg) = ai^E^’m(s)- 

It can be seen that ayf^^’m(s) = 0 for i 7 ^ j, gy^^$m(s) = 0 for i < j, A: < with ij ^ kh, oyf^^m(s) = 
0 for j < j, k < h, and gy^.y^^ ^ 77 i(s) = 0 for i < j, A; < /i with ij ^ kh. 

Using the previous computations, the operator m applied to $m(s) gives 


q 2 d. q2 

E ,kh p. ^ ^ -^kk,kk q 2 ^'A7i(-Es) 

^^ijZkh rinr'. . 


{i,j){k,h) 


k=l ^^kk 


Q2 q2 

^ Akh,kh ( ) $m(3fs). 




9xlh dylf^ 


and now from (1161) and (d3 we get 


( 22 ) 


^ ^ q2 ^ ^ ^ 

^ ^ -^kk,kk ^ ^ ^ -^kk^kk ^ ^ 

k—1 


I'^km'f^kj r 


k—1 


(23) 


and 


/ g2 q2 

E ) $m(3fs) 

l<k<h<d \^^kh ClVkh 


(24) 


— 2 ^ ^ -^kh,kh ^ ^ 


l<k<.h<d 


j^m 


^ I'^km'^hj '^hm'^kjl ^ I'^km'^hj '^hm'^kjl 


Am Ad 


Am Ad 


2 ^ ^ -^kh,kh ^ ^ 2 


\Ukm(s)Uhj{s)f + \Uhm{s)Ukj{s)f 

Am ( 5 ) Aj (s) 


Finally, dUD follows from 


Remark 16 (a) Theorem \15\ is the analog of the Dyson-Brownian motion {Ip for a Hermitian Levy process. 
In this case, X{t) has a simple spectrum almost surely for each t > 0. The problem of noncolliding eigenvalues 
for each A > 0 almost surely will be considered elsewhere. 
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(b) However, we observe that the repulsion term between each pair of eigenvalues l/{Xi{s) — Xj{s)), i ^ j, 
appears in m and 117]) when there is a Gaussian component (i.e., A^Q). 

(c) The corresponding expression of Proposition^ can be obtained for the dynamics of eigenvalues in 

terms of the local martingale and bounded variation components of the eigenvalues. Let {Xt : t > 0} be a 
Hermitian Levy process with triplet [A,Vt'^) satisfying the same conditions as in Theorem 1751 For each 
m = 1,2, it holds the decomposition Xm{Xt) = where is a martingale and is a 

process of bounded variation such that 


= Xm {Xo) + tr ^ DXm {Xs-)BA{dsfj 


/(0,t]xH« 


[Xm [Xs- +y) - Xm Jx {ds, dy) 


and 


yr = \l 


{i,j){k,h) 


ij ^kh 


Xm{Xs)ds + tT^J '^DXm{Xs)ds'^ 


+ / [Xm (Xs- +y)- Xm (Xs-) - tr {DXm (Xg-) y) l{||y||<i}] dsv{dy). 

7(0,t]xH« 

d) Hermitian Levy processes with rank one jumps and with their distributions invariant under unitary 
conjugations appear in the context of a connection between classical and free one-dimensional infinitely 
divisible distributions /73|/. gg/, m/, \23^ . Let us consider an example with the triplet [ 0 ^ 142 , i/, ijl. From 
(c), for each m = 1,..., d, the eigenvalue Xm has the decomposition Xm{t) = , where 


= Xm {Xo) + tr ^ DXm {Xs-) BA{dsfj 

+ / [Am {Xs- +u)- Xm {Jx {ds, du) - dsv{du)) 

7(0,tlxHi 


and 


vr = 2am 


0 j^m (^s) 


ds + '0tr^y DXm{Xs)ds 


+ / [Xm{Xs-+u)-Xm{Xs-)-tT{DXm{Xs-)u)l^\\u\\<i}]dsr{du). 

J (0,t] XH;^ 


In particular, when the process is symmetric, the Levy measure is symmetric and -0 = 0, and hence 


vr = 2r 


1 


0 j^m (Xs) 


ds 


+ [Xm{Xs-+u) - Xm{Xs-)]dsiy{du). 

7(0,i]xHi 
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